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Abstract 

We investigate solutions to a nonlinear integral equation which has a central 
Q\ \ role in implementing the non-Abelian Gauss's Law and in constructing gauge- 

invariant quark and gluon fields. Here we concern ourselves with solutions to 
"q 1 this same equation that are not operator-valued, but are functions of spatial 

variables and carry spatial and SU(2) indices. We obtain an expression for 
Q\ , the gauge-invariant gauge field in two-color QCD, define an index that we will 

refer to as the "winding number" that characterizes it, and show that this 
| winding number is invariant to a small gauge transformation of the gauge 

field on which our construction of the gauge-invariant gauge field is based. 
We discuss the role of this gauge field in determining the winding number of 
the gauge-invariant gauge field. We also show that when the winding number 
of the gauge field is an integer ly^O, the gauge-invariant gauge field manifests 



o 

Q\ ' winding numbers that are not integers, and are half-integers only when I = 0. 

ON " 



I. INTRODUCTION 

In earlier work, in which we discussed QCD in the temporal (AJ = 0) gauge 



X 



| we 

implemented the non-Abelian Gauss's law and constructed quark and gluon operator-valued 
fields that are invariant to small non-Abelian gauge transformations. An essential element 
in that program was a nonlocal operator, A](r), which we called the "resolvent field" in 
later publications. [|2]|| The resolvent field is a functional of the gauge field A](r), and it 
has a pivotal role in our work. It first appears, in combination with U] (r) — the canonical 
momentum for the gauge field and the negative chromoelectric field — in the operator, 



» = ||«P W *^(r W (,)j||, (1.1) 

where the double-bar bracket denotes a normal ordering in which all the gauge fields are 
placed to the left of all canonical momenta, so that when the exponential is expanded, 
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all powers of A i (r) appear to the left of any of the n/(r). ^ converts any state \(j>) that 
obeys the Abelian Gauss's law djl" (r)|0) = into a state ^\4>) that obeys the "pure-glue" 
non-Abelian Gauss's law [<%nf(r) + gf a ^A? (r)ir/(r)]^|0) = . Q 

The resolvent field reappears, in the form 

3^(r) = %Aj{v) , (1.2) 
in the expression for the gauge-invariant quark field 

^Gi(r)=Vc(r)V(r) (1-3) 

where 

V c (r) = exp(-i^(r)f ) exp (-^ a (r)f ) (1.4) 

and 

AT°(r) = [^(r)]. (1.5) 
The resolvent field also appears in the expression for the gauge-invariant gluon field 

[Aln(r) f ] = V c (r) [AJ(r) f ] V c ~\v) + J lfc(r) ^(r) , (1.6) 
which can also be expressed as 

A]{v). (1.7) 

In subsequent work, [^-^] we used these gauge-invariant fields to derive an expression 
for a nonlocal interaction that couples gauge-invariant quark color-charge densities to each 
other and to gluons through "chains", in which each link consists of an integral over the 
gauge-invariant gluon field. SU(3) or SU(2) structure constants connect the links to form 
the chains. As was shown in Ref. |]J, the quantities that appear in this formulation - 
the gauge-invariant quark color-charge density j'q G | = ^gi^V'gi and the gauge-invariant 
gauge field Aq U {v) — commute with the generator of infinitesimal gauge transformations 
and therefore are completely invariant to small gauge transformations. In the work presented 
here, we will focus attention on the resolvent field in two-color QCD. We will examine its 
structure in configuration space and its role in determining the topological features of the 
gauge-invariant gluon field. 



1 We use non-relativistic notation in this work, in which subscripted indices designate contravari- 
ant components of contravariant quantities, such as Aj , and covariant components of covariant 
quantities, such as d{. This convention is used extensively throughout this work. 
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II. CONFIGURATIONS OF THE RESOLVENT FIELD 



In Ref. we presented a nonlinear integral equation for the resolvent field. Although, in 
Ref. our concern was with the series representation of the operator- valued resolvent field, 
it is also relevant to ask whether the integral equation for the resolvent field has solutions 
other than that series. It is this question that we now consider. The most general form 
of the integral equation for the resolvent field, applicable to any SU(N) gauge theory, is 
not useful for this investigation; but its SU(2) form, given in Ref. [Ill], is very suitable for 
our purposes. In the SU(2) case, the integral equation for the resolvent field A](r) can be 
expressed in the following form: 



where 



Al(r)=AJ(r) x + Al(r)y, (2.1) 



_ ge «M Q. X P (T) 

_ g 2 e aPb ebfil x »( T j x a( T ) 1-oob(A0 



and 



Al(r)y = ge^y*(r) (A^(r) + (* y - ^)A%v) ) ^fQ 



+5^ 1 f(r)a i F(r)^ 



7T 



+g h^ b e b ^ y»(r) y«(r) ^(r) - ^] . (2.3) 

We observe that the right-hand side of Eq. ([2.2|) is independent of A](r) and depends only 
on the gauge field A](r) and on trigonometric functions of Af, where 

A" = Wx a {r) X a {v)Y ■ (2.4) 

Aj{r)x therefore is not responsible for any of the nonlinearity of the nonlinear integral equa- 
tion shown in Eq.( |2.1| ), but represents an inhomogeneity in that equation. The right-hand 
side of Eq.( [2.3|) , however, contains the resolvent field AJ(r) and AT as well as trigonometric 
functions of A*, where 

Jf= [g 2 y^{r)y^{v)f . (2.5) 

Aj{r)y therefore is a nonlinear functional of the resolvent field A](r), and of 3^ 7 (i"), whose 
dependence on A](r) is given in Eq. ( |1.2| ). It is the structure of A](r)y that makes Eq. 
a nonlinear integral equation. 
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Eqs.( |1.3|) and ( |1.6|) are identical in form to the equations describing the change in a spinor 
field and a gauge field, respectively, produced by a finite non-Abelian gauge transformation, 
with the important exception that the transforming c-number gauge function in a gauge 
transformation is replaced by the operator- valued gy a (r) and gX a {r) in Ve(r). The SU(2) 
version of Vc(r), given by 

[Vc(r)] sm = exp ( -igy^{vYi ) exp (-i^(r)^) , (2.6) 

can be interpreted as a sequence of two rotations represented in the fundamental represen- 
tation of SU(2). In one of these rotations, gX a (r)/Af designates a direction about which the 
rotation is made, and M the magnitude of the locally varying angle of rotation; a similar 
correspondence applies to gy (r) and 77 in the other rotation. In this SU(2) case, we can 
also express [Vc(r)] 5[/ ^ in the form 

[Mr)W (2) = exp ( -igZ^{vYi ) (2.7) 

where the composition laws for rotations enable us to express Z a (r) as a function of X a (r) 
and y a (r), and both in terms of equivalent sets of locally varying Euler angles, much as was 
done by Christ and Lee. Because the algebra of SU(N) generators is closed, the SU(3) 
form of Vc(r) can also be expressed in the form Vc(r) = exp igZ a (r)^Y^j , where the 

Baker-Hausdorff-Campbell theorem enables us to equate Z a to a complicated expansion in 
X a and y a . The familiar group composition laws for rotations are a special case of the 
Baker-Hausdorff-Campbell theorem applied to the SU(2) case. Eqs.( [2.2|) and ( |2.3|) illustrate 
this similarity in form of Eqs.( |1.3j ) and ( |1.6| ) to the equations that implement finite gauge 
transformations. Eq.( [2.2|) has the same form as a finite gauge transformation of a gauge field 
in the adjoint representation of SU(2), with gX a (r) serving as the function by which the field 
is gauge-transformed; the same applies to (|2.3|) with gy a (r) in place of gX a (r). Such non- 
Abelian gauge configurations have well-documented important geometric and topological 
implications for gauge theories. [|7|-|T0H The similarity in form of the equations that establish 
gauge-invariant fields to those that effect gauge transformations accounts for the appearance 
in our work of expressions that resemble the inverse of the Faddeev Popov operator di(S ac di + 
gf cba A b ), which arises when transformations from the temporal to the Coulomb gauge are 
carried out by treating them as transformations from Cartesian to curvilinear field variables. 
P JTT|JT2]] In our work, di(S ac di + gf cba A b Q H ) appears in place of the Faddeev-Popov operator, 



and it appears in a representation in which the gauge field A b QU and the quark color charge 
density jg = gip\)^ /2)ip are gauge invariant and in which the temporal gauge condition 
continues to apply, even though the gauge-invariant gauge field (but not the gauge field) is 
transverse. j|,|5| 

In spite of this similarity in form of finite gauge transformations to Eqs.( |1.3| ) and ( |1.6| ), the 
significance of these two equations is very different from that of a finite gauge transformation. 
In a gauge transformation, the gauge field is transformed by an arbitrary c-number function 



that bears no relationship to the original gauge field. But in the case of Eqs.(L3) and 
(|OD, X a and y* (as well as Af and AT) are themselves functionals of the gauge field, and 
therefore are subject to the same gauge transformations that affect it. In fact, under the 
infinitesimal gauge transformation SA](r) = <9i<5u; 7 (r) + ge^^A 1 * {y)5ujP {r) , Vc(r) transforms 
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so that 5Vc(r) = —igVc(r)5u' y (r)^-. [[L3[] Instead of representing a gauge transformation, 



Eqs. (|1-3|) and flL6| ) establish the operator- valued nonlocal quantities ipG\{ r ) an d Aq H (t), 
which are entirely invariant to any further transformations effected by the generator of 
gauge transformations 



exp i J dr ^nj(r) + gpf>«A? (r)Uf (r) + #t( r )^( r ) 



w 7 ( 



(2.8) 



The integral equation described in Eqs. (|2.1| )- (|2.3| ) can be viewed in two distinct ways: On 
the one hand, as was shown in Ref. [Q, Eqs. (|2.1|) -(|2~3|) can be used to generate an operator- 
valued series representation of the resolvent field A](r), which is instrumental in establishing 
gauge-invariant quark and gluon fields and quantum states that implement the non-Abelian 
Gauss's law. But, because all the operator-valued fields that appear in it commute with 
each other (the only operator that would not commute with the quantities in Eqs.( [2.1|) - 
( |2.3D - - II"(r), the momentum conjugate to Af(r) - - does not appear in any of these 
equations), Eq. (|2.1|) also can be treated as a non-linear integral equation whose solutions 
can be investigated by more-or-less standard procedures. In this second context, we can 
view the resolvent field A](r) and the gauge field A](r) as functions of spatial variables that 
obey the same nonlinear integral equation as the corresponding operator-valued quantities. 
We can, furthermore, make an ansatz about the functional dependence of A](r) and A](r) 
on spatial variables that is physically reasonable, and that facilitates the solution of the 
nonlinear equation given in Eqs.(2.1)-( |2.3| ). The functional forms of A](r) and AJ(r) will 
represent number-valued realizations of these operator-valued fields, and will enable us to 
investigate their spatial configurations and their topological features. 

We will carry out this program by assuming that A](r) and ^47( r ) are functions of spatial 
variables, that they are second-rank tensors in the combined spatial and SU(2) indices i and 
7 respectively, and that, except in so far as the forms of *47( r ) and A](r) must reflect this 
second-rank tensor structure, they are isotropic functions of the position. The most general 
form of A] (r) will therefore be 

AJ(v) = S^Ur) + rj ^ L Ur)+e i ^Ur) (2.9) 

where ^(r), an d £c( r ) represent as yet unspecified isotropic functions of r. It is con- 

venient to separate the resolvent field A](r) into longitudinal and transverse parts, A] L (r) 
and A] T (r) respectively. We then use the representations 

7J L (r) =5 i7 $(r) + ^$ / (r) where $'(r) = ^ (2.10) 



and 



AJ T (r) = 6 tl <p A {r) + ^ Mr) + ^J- Mr) (2.11) 



nn2 " ' ' iy 



where the transversality of A] (r) requires that 



d(r > ) +r 2 ^ = 0. (2.12) 
dr dr 
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The transversality of A] T (r) does not impose any conditions on (pc- On the basis of these 
considerations, we can write the resolvent field as 



A](r) = 5^ {ip A + $) + ^ (<p B + r$') + e i7 „^ , 



and y a (r) becomes 



y a (r) = r Q $(r) 



(2.13) 



(2.14) 



The same analysis can be applied to the gauge field A](r), with the result that we obtain 
the wholly analogous equations 



and 



A] T (r) = 5^ T A (r) + ^ Mr) + ^J- T c (r) 



(2.15) 
(2.16) 



with 

d(r 2 T B ) | r2 dT A 
dr dr 

The complete gauge field can be represented as 



gr 



A/" + 



r 4 r 7 



cfAT 



r dr 



+ Ta H — Tb + ej 7 „ — Tq , 



(2.17) 



(2.18) 



where we have used Eqs. (|2.4|) and (|2.15|) to set grS = 77. 



We can use the representations given in Eqs.( p.lO| ), ( |2.11|) , (|2.15|) , and ( |2.16|) to express 
Eq.(O) as 



<A]y = \ e hn — [p> A + Ta] - ( S i7 - ^-f- ) [ip c + T c ] \ sin AT 



+ — e i7 n — ( 1 — cosA/* 



gr r 



n r 7 



+ L( 5 . -Till 

gr V r z 



[(fi A + Ta] + e iyn — [cpc + T C ] \ (l- cos AT 



(77- sin 77 



(2.19) 



where we have used Eqs.( j2.5p and ( p. 14 ) to set gr$ = 77. In a similar way, we can express 
Eq.(p|) as 



A 1 — 



TV 

r 



«7n 

gr r 



77 ' 
— + T A 
5-r 



1 - cos A/") 



"*7 r 2 



TcUinA/" 



7 r 2 



A/" ^ 

5-r 



5-r 



'7 r 2 



+ e i7n yT c | (1 - cos A/") 
(AT -sin AT) . 



(2.20) 
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It is manifest from Eqs. (|2.19|) and (|2.20|) that r 1 A 1 i {r)-^ = and r 1 AJ(r)y = 0, so that 



y A](r) = too. When we use that fact in Eq. (|2.13|) , we find that 

AT' 



$ + r$' + Lp A + ip B = or, equivalently, \-(p A + ip B = 0, (2.21) 

9 



enabling us to represent A](r) as 



Al(r) 



Si 



7 T^T 



K gr 



+ ¥a 



■<pc- 



(2.22) 



Substitution of Eqs. (|2.19| ), (|2.20| ), and ( |2.22|) into Eq. (|2.1|) enables us to express <p A and ipc 
as functions of Af, Af, T A , and T c , so that we obtain 



= — 

gr 



Afcos(AT + A/") - sin(AT + Af)] + T A [cos(AT + M) - ll - T c sin(AT + Af) (2.23) 



and 



gr 



A/"sin(AT + Af) + cos(AT + Af)-l]+T c [cos(AT + AT) - ll + T A sin(AT + Af ) . 



(2.24) 



We can solve Eqs.( |2.12j ) and ( |2.21| ) simultaneously to eliminate (fB, and obtain a nonlin- 
ear differential equation that is equivalent to the nonlinear integral equation described in 
Eqs.( |2.1| )-( |2.3j ) modulo the previously- mentioned assumptions about the forms of the re- 
solvent field and the gauge field. It is most convenient to use the dimensionless variable 
u = ln(r/r ), where r is an arbitrary constant length, as the independent variable in ex- 
pressing this differential equation, which then is 



d 2 Af dAf 
+ 



du 2 du 

+ 2#r exp(u) \j A 



A/"cos(A7 + Af) - sin(AT + Af) 



cos(AT + Af) - l] - T c sin(AT + Af)} = 0. 



(2.25) 



Eq. ( |2.25|) is a nonlinear differential equation in the variable Af, which determines the resol- 
vent field. Af, Xa, and Tq will be taken to be known functions in this equation. One question 
of immediate interest is whether there are nonvanishing solutions for the case that Af, T A , 
and T c all vanish, so that the gauge field A] is identically zero. In that case, Eq.( |2.25D 
reduces to 



d 2 Af dAf n . ^ 
+ — 2sm(AT) 



du 2 



du 



0. 



(2.26) 



We observe that Eq. (|2.26|) is the equation for a damped pendulum in which small-amplitude 
oscillations have not been assumed, in which u corresponds to the time, and Af — 9 + ti, 
where 8 = when the pendulum is in its stable equilibrium position. Moreover, Eq.( [2.26 ) 
corresponds to the equation for the pendulum in which mass, length, local acceleration due 
to gravity, and damping constant have been set = 1, so that no adjustable parameters remain 
in the equation. 
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In discussing Eqs.( |2.25| ) and ( |2.26|) , we will assume that T4, Tc, and Af — 2£n (where £ is 
an integer), vanish faster than 1/r as r— >oo. Moreover, we also assume that the gauge field 
is bounded everywhere, so that rT^— >0 as r— >0, and that the same applies to rT c , and rAf. 
We will make these conditions more precise in the Appendix. 



The fact that Eq.( [2.26 ) describes a damped pendulum moving without a driving force 
makes it obvious that it has nonvanishing solutions. In the limit u^oo, the damped pen- 
dulum must come to rest in a position of static equilibrium, so that, as u^oo, A/*— >rrm, 
where m is an arbitrary integer. Physical considerations lead us to expect that the damped 
pendulum will come to rest in a position of stable equilibrium, for which A/*— >{2m + \)it. 
But in the application of Eq.( 2.26| ) of interest to us in the current context, all possible so- 
lutions — those that terminate at saddle points as well as those that terminate at stable 
equilibrium positions — must be considered as long as they are bounded functions in the 
interval — 00 < u < 00. Eq. ( |2.25| ) is more complicated than Eq. ( j2.26| ); but as a result of the 
conditions on T4, and T c , and on Af for the £ = case, discussed above and made precise 
in the Appendix, Eq. ( |2.25| ) approaches Eq.( |2.26D at the origin in configuration space, i. e. 
in the limit u—>- — 00, as well as in the limit u^oo. It is therefore reasonable to expect that 
A/*— >rmr as u — >oo even when it is governed by Eq.( |2.25D . That this expectation is fulfilled, 
is shown in the Appendix. When £7^0, A/* still approaches a definite limit as but that 

limit is no longer an integer multiple of 7r. This case is discussed in Section OIL and the 
mathematical foundations for this discussion are given in the Appendix. 



It is possible to use the results obtained in this section, in particular Eqs.( |2.H] ), ( [2.16D , 
and flpD, to represent the gauge-invariant gauge field as 



.4 



Gl , 



1 

gr I r 



cos(AT + A0 - 1 +A/"sin(A7 + A0l + [6. 



'1 7 



x 



Afcos(AT +Af) - sin(AT + Af) 

TnT, 



+r A 
+r c {e 



Oj 7 



r r . 



1' 7 



cos (AT + Af ) 



r i r 1 dAf 
r dr 
r 



c y ,- n — cos(AT + Af) 



1 7 



sin(AT + Af) 
sin (AT + Af) 



(T A + T B ) 



When we set the entire gauge field A] = 0, we obtain 



[AIM 



(AJ=0) 



gr 



-1711 



cos 



at; 



Si 



17 



sin. 



Af 



rjr 7 dfN 
r dr 



(2.27) 
(2.28) 

7 ■ wi j 

As a consequence of Eq.( [2.25| ), diAQ H (v) = 0, and similarly, di [^4Gu( r )](A 7 =o) = follows 
from Eq. (|2.26|) , so that the consistency of Eqs.( [2.27| ) and (|2.28| ) with the transversality of 
the gauge-invariant gauge field is confirmed. 

In order to investigate the changes in AT that can ensue from gauge transformations, we 
turn our attention to the changes SAf, ST A , STb, and STc in the corresponding quantities 
Af, T4, Tb, and Tc due to the infinitesimal SU(2) gauge transformation 



SA](r) = «W(r) + ge ial3 A?(r)SujP(r) 



(2.29) 



S 



We extend the ansatz that led to Eqs. fl2.15D and ( 2.16|) to choose an expression for 5uj 1 {y) 
that is isotropic except in so far as it transforms as a vector in the adjoint representation of 
SU(2). We therefore represent Suj 1 (r) as 



5^{r) = ^A(r) 
r 



(2.30) 



where A(r) is an infinitesimal, isotropic quantity. We substitute SAf, ST A , 5T B , and 5Tq for 
Af, T A , Tb, and Tq respectively, in Eq. (|2.18|) , to represent 5AJ(r), the infinitesimal change 
in A](r) due to a gauge transformation, and obtain 



5i 



17 



Ti r 7 \ 5 AT Ti r 7 d (SAf) 
gr gr 2 dr 



ty . ty ry 

ST B (r) + e lin ^ 5T c (r) . 



(2.31) 



We then substitute Eq.(|3T]) in the left-hand side of Eq. (fO|), and Eq.(gT§) in the right- 
hand side to obtain 



5Af(r) 
gr 



ST A (r) 



A(r) 



gA{r)T c {r), 



(2.32) 



d (5N{r) 



dr \ gr 



+ ST B (r 



A(r) 



d_ 

dr 



A( r ) - gA(r)T c (r) 



and 



ST c (r) = -gA{r) 



'JV(r) 



gr 



+ T A (r) 



We combine Eqs.([2TT7|), and (|2733|) to obtain 

d 2 „ . , , _ d 



r 2 — 50(r) + 2r— <50(r) - 2<50(r) = -2#rA(r)T c (r) 



fir 



where we have set 



dr 



5Af(r) 
9 



A(r) = 6<i>(r). 



With a little more algebra we can also obtain 

r d 2 



ST A (r) 



d 

dr 



and 



r d 2 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



To investigate the changes that gauge transformations induce in Af, we examine Eqs. ( ^.23| ) 
and ( |2.24j ) to find relationships between the changes that infinitesimal gauge transformations 
produce in the quantities that appear in these two equations. We find that 
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5^ A = — [577cos(77 + 77) - [A/"sin(AT + M) + cos(AT + M)] (577 + 577)} 



gr 
5T A 



cos(AT + 77)-l]-T A sin(AT + A/") (577 + 577) 



- ST C sin(AT + 77) -T cos(AT + 77) (577 + <5A/") 



(2.39) 



and 



5 if 



<= |<5jVsin(A7 + A/") + [a/"cos(A7 + A/) - sin(AT + A/")] (MT + &/V)} 

+ 5T C [cos(AT + AT) - l] - T c sin(AT + 77) (577 + &/V) 
+ 57 A sin(AT + A/") + T A cos(AT + A/") (5AT + 577) . 



^2.40) 



As shown in Eq. (JL7| ), changes produced by gauge transformations in the transverse compo- 
nents of A] and A s cancel, so that 



5ipA + 5T A = 0, 5(f B + 5T B = 0, and 5(p c + 5T C = 0. 
We substitute Eq.(ggl]) in Eqs. (ggg) and (fODD , set 

X = T C + — , F = T yl + — , = A7 + AT, and £ = 5AT + <W - #A 



5-r 



and find that 



(2.41) 



(2.42) 



(ATcosfl + Fsinfl) £ = 



(2.43) 



and 



(ATsin# - Ycos9) £ = 



^2.44) 



Since X and Y are determined by the initial configuration of the gauge field A], which is 
arbitrary, Eqs. fl2.43]) and ( |2.44j ) must be satisfied without requiring that X and Y vanish. 
It therefore follows that £ = 0, which, with Eq.( 2.36| ), leads to 



577 



(2.45) 



and, with Eq. (|2.35| ), to the following differential equation for the variation of 77 produced 
by an infinitesimal gauge transformation: 



d 2 

r 2 —577(r) 
dr z 



2r^-577(r) - 2577(r) 
dr 



2g 2 A(r)rT c (r) 



(2.46) 



Further confirmation of the consistency of these results is obtained by examining the varia- 
tion of Aq H (t) due to an infinitesimal gauge transformation of the quantities that constitute 
it. This variation is given by 
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1 

gr 



-1 7 n 



5M sin(AT + M) + (ATcos (AT + AO - sin(AT + M)) 5 (AT + M) 



+ 



'17 



<5A/"cos(AT + AT) - (ATsin (AT + AT) + cos(AT + AT)) 5 (AT + A/" 



r;r 7 d(5AT) 



(/?• 



r;r 7 



1 7 



-j 7 n 



cos(AT + A/") + e 47n -sin(AT + AO \ + ^ (5T A + 5T B 



-cos(AT + AT) 



* 7 



r. 



sin (AT + A/") 



sin(AT + M) - e iin — cos(AT + AT) U (AT + M 



-To {e l7n ^sin(AT + A/") + f<^ 7 — ^) cos(AT + AT) } 5 (AT + Af) . 



(2.47) 



When Eqs.(P35D, ( P^op , ( gg7| ), ( ggg ), and ( P^BQ are used to replace 5A0 <5AT, 5T A , 5T S , 
and 5T C in Eq.( [2.47| ), we find that <5Ag U (r) = 0, confirming the consistency of the procedure 
that leads to Eq.( ggE| ). 

The effect of small gauge transformations on the limits of 577 as u— > ± oo can best be ex- 
plored by expressing Eq. Q2.46Q in terms of the independent variable u. As is discussed in the 
Appendix, the expression rTc appearing on the right-hand-side of that equation is subject 
to the conditions lim^oo |r exp^T^u) | < K exp(-au) and lim<^_ 00 |r exp(u)T c (u) \ < 
K' exp(cra), where K and K' are constants. Constraints can be placed on A, but the only 
one necessary to complete this argument is that A is bounded in the interval (— oo, oo). 
These circumstances allow us to conclude that Eq.( 2.46 ) can be represented, in the limits 
u^oo and u— * — oo, in the form 



du z du 



< C exp(— ot\u\) 



(2.48) 



where C is a constant. We can therefore infer that 

lim \57T{u)\ < Aexp(au) + £>exp(u) and lim |5AT(m)| < v4'exp(— era) + B'exp(— 2u) 

(2.49) 

for some constants A, B, A', and B' . We observe that the limiting values of 5J\T(u), as 
tt— >±oo, vanish, and that therefore the corresponding limiting values of A/*(u), as u— >oo and 
as u — > — oo, are invariant to small gauge transformations — i. e. transformations implemented 
by diU](r) + gt^^A^ (r)II"(r) + gip^ (r) 1 ^-^^) acting as the generator of infinitesimal gauge 
transformations. But AT(u), for arbitrary values of u, does not share that invariance. These 
arguments will be used, in the next section, to show that the winding number of the gauge- 
invariant gauge field remains unaffected by small gauge transformations. 



III. TOPOLOGY AND THE IMPLEMENTATION OF GAUGE INVARIANCE 

Our considerations in this section are based on the fact, discussed in Section 0, that 
the transformation that converts the quark field ip(r) and the gauge field A](r) into their 
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respective gauge-invariant forms — as shown in Eqs.( |1.3|) and (|1.6|) respectively — has the 
formal structure of a gauge transformation, although it differs from a gauge transformation 
in a very important respect: The argument of the operator- valued Vc(r) that effects the 
transformation to gauge-invariant forms is not a function by which ip(r) and A](r) are gauge- 
transformed, but is itself operator- valued and subject to gauge transformations. Eqs.(Ol) 



and (|1.6|) therefore actually represent transformations of gauge- dependent fields to forms that 
are invariant to further gauge transformations. Nevertheless, when we represent Vc(r) as a 
number-valued realization of the operator-valued quantity, we can use its formal structure 
to investigate the topological features of the gauge-invariant gauge fields. 

With the ansatz that led to the expression for the resolvent field given in Eq. ( |2.22j ), we 
find that Vc(r) is given by 

(77 + Af) 

V c {r) = exp I -if n rj >- I . (3.1) 



and that the expression for AQ H (r) given in Eq.( |2.27 ) can be written as 



w\Al H (v) = -igV c {v) T -A1{v)Vs l {v) + VcmVc'ir), (3.2) 



where A](r) is given in Eqs.( p.l8|) , and where we can define 



A<(r) = -ig- [Al H {v)\ v = ^(r^V^r) (3.3) 



so that [Ag H (r)]y is the part of Ag U (r) given by 



[AIM] 



v ~ (cos(A7 + A/) - 1) - (* 7 - ^) sm(A7 + A/) - ^ ^±^j 

(3.4) 



i' 7 



Eq.( |3.3|) shows that [^guWIv nas the formal structure of the "pure gauge" part of a gauge 
field, although that is not what it is. [Aq\ j(r)] v is a part of the gauge-invariant field ^^(r), 
and any change in [^4Qi( r )]y due to an infinitesimal gauge transformation must be offset by 
compensating gauge transformations to Vc^^AJ^V^ 1 ^). We can, nevertheless, make 
use of the formal structure of Aj(r) represented in Eq.(|3.3|) to define 

Q = -(24^)- 1 ^/*lV[A ( (r)A J -(r)A t (r)]; (3.5) 

and, following Ref. 0, we can use Eqs.fl3.lp and ( |3.3p to express Q as 



_ 1 fdv , 2 ( J\F + Af \ d (aT + A/") 

Q = iW [—2- ) dr ' (3 ' 6) 



and integrate this expression to obtain 
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Q = — Urn (AT(r) + Af(r) - sin (AT(r) + Af(r)) } 
{^(0) + AT(0) - sin (AT(0) + Af(0)) } . 



(3.7) 



We will refer to Q as the winding number of the gauge-invariant gauge field, even though 
its values are not restricted to integers. Unlike the winding numbers of the gauge field, 
which are integer-valued, |7HT0| Q is determined by the properties of 77, which is governed 
by Eq.( |2.25| ) and by the fact that it must be bounded in the entire interval 0<r < oo. In 
discussing the properties of Af, we define the "pure gauge" gauge field 



y 47( r )pg 



1 

gr 



-1711 



[cos(JV) - 1] - [Si 



r i r 1 



sinAf-^^ 

r dr 



which we obtain from 
AKr 



pg 



-Tr 



exp ( —iJ7 



~2r~ 



di exp ( iAf 



~2r 



(3.8) 



(3.9) 



We further observe that, in Section |TJ and in the Appendix, we have assumed that Af = 
at r = 0, and that lim^oo Af = 2£tt, so that we consider gauge fields A](r) pg whose winding 
number is the integer £, and therefore that 



Q = — <^ lim 77 (r) - sin (77 (r) 



2n£ - 77(0) + sin (77(0) 



As is shown in the Appendix, when £ = and 77 is normalized so that 77 = when r 
lim^oo 77 = mix. For the £ = case, Q is therefore given by 



(3.10) 
0, 



Q = — ( lim Mr) - 77(0) 

2tt y-»°o 



— lim 77(r) 

2tt r -*°° 



m 
~2 



(3.11) 



where m is an integer. We will now discuss the implications of Eq. ([2.25 ) for the values of Q 
when £^0. 

When we consider gauge fields with winding numbers £^0, for which lim^oo Af(u) = 2tt£, 
but limu^oo exp(u)T^(u), and lim^oo exp(u)7c(u) still vanish along with liniu^.oo Af(u) , 
lim^-oo exp(v 
77(u) satisfies 

d 2 77 d77 



lim u ^_ OQ exp (u)T^(u) , and lim u _^„ 00 exp(u)Tc(u), then we observe that, in the limit u— >oc, 



du 2 



du 



27r£cos(AT + 2n£) - sin(AT + 2 



Tit 







(3.12) 



and 77(u) approaches a limiting value lim^oo 77 (u) , which exists, but is not an integer mul- 
tiple of 71. 77 (u) vanishes, as before, at the saddle point for which u— * — 00, but lim n _^ 00 77 (u) 
now is one of the denumerably infinite solutions of the transcendental equation 



tan < lim 77 ( 



u, 



2tt£. 



(3.13) 



If we adopt the convention that the inverse tangent is defined so that, for £ = tan 1 (rj), 
— I < £ < |, then the winding numbers for the gauge-invariant gauge field are given by 



13 



Q = i- {tan- 1 (2 7 r£) + mvr + 2^ f 1 - } (3-14) 



for some integer m. 



Eq.( |3.14 ) establishes a winding number for the gauge- invariant gauge field that is a 



function of two integer-valued variables, I and m, but that is not itself integer-valued. 
The integer £ defines the homotopy class of the gauge field described by Eqs. (|3.8|) and 
(|OD. Gauge- invariant gauge fields can be categorized by the value of t that describes the 
homotopy class of the gauge fields to which they are linked by Eq. (|2.25| ), and by the value 
of m that defines the sheet on which the inverse tangent function that corresponds to the 
lim^oo AT (it) is located. The range of values of the integer m, for which mir corresponds to 
limits that Af(u) can actually attain, depends on the solutions of Eq. ( |2.25| ) that given sets of 
functions Af(u), Ta(u) and T c {u) can support. Examples of solutions, and their implications 
for winding numbers of the gauge-invariant gauge field, are given in Section [TV| 

The relationship between Q, I and m reflects the fact that, when lim u ^_ 00 M{u) = 0, the 
lim^oo M{u) is the only part of the gauge field that can affect its homotopy class; however, 
the functional dependence of Af(u), Ta(u) and T c {u) for the entire range of values of u affects 
the lim^oo AT(-u), and therefore Q, as shown in Eqs.( p.25|) and (|3.10|) . A number of different 



values of m can therefore be compatible with the same £, and, in fact, with the identical 
functions A/"(it), T^iu) and T c {u). The mathematical foundations for these observations are 
discussed in Section [V] and in the Appendix. 



IV. ILLUSTRATIVE NUMERICAL EXAMPLES 



In this section, we present illustrative examples of numerical integrations of Eq.( f2.25|) 



that simulate solutions that are bounded for all real values of the independent variable.^] In 
carrying out these numerical integrations, we have selected functional forms for Ta, Tq, and 
J\f that are consistent with assumptions ( |A1| ) to (|A3| ) and that include examples of M that 



belong to the trivial homotopy class with i = 0, as well as others in which lim^oo M = 27i£ 
with £y^0. The technical considerations on which these simulations are based, and the 
individual examples, constitute the remainder of this section 

In Theorem in the Appendix, we have established the existence of a solution AT sat- 
isfying equation (|A4|) and boundary conditions ( |A5|) , provided assumptions ( |A1|) to ( |A3|) 
are satisfied. (It is noted that w and y in the Appendix represent A^ and AT, respectively.) 
Such assumptions essentially say that e*7^ and e l Tc behave like for large |t|, and that Af 
behaves like as t — >• — oo, and like llix for some integer £, as t — > oo. 



2t 



In this section, and in the Appendix, Eq. fl2.25 ) will sometimes be written in the form of Eq.(A4), 



and the independent variable u will sometimes be written as t. 
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Remark 3 in Theorem [1] asserts that there are an infinite number of solutions to Eq. 



But the question arises: Will numerical solutions of Eq.(S4), that are bounded in the entire 
interval (—00, 00) and normalized so that lim^.oo AT = 0, always settle down to the same 
two values of m as in the case for the damped pendulum? (It is well known, in the case of 



Eq. ( |2.26| ) which describes a damped pendulum, that there are exactly two classes of solutions 



- one for m—1 and another for m = —1.) We will carry out some numerical experiments 
to address this question, and document some of our results in this section. 

Since our numerical computations will be done on a finite interval [-M, M] for some 
large value M > 0, we will first construct the appropriate boundary condition at £ = —M. 
With the solution AT going to zero as t — > —00, one can impose the boundary condition 
Af(—M) = 0. However, a better boundary condition can be constructed, which leads to 
reduced numerical errors even when we use only a moderate value of M. This reduced 
interval size for [— M, M] leads to more efficient and more accurate calculations. Such 
attributes of numerical infinity are well known (see Chapter 4, Ref. P~j). 

Since both AT and M tend to zero as t approaches —00, we employ a Taylor's expansion 
to extract the leading order behavior of the terms in equation ([S4]) in that regime. We keep 
only the leading order terms in Af and Af, taking into account that while both are small 
when t— > — 00, the relative magnitudes of Af and Af are not known. In this way, we find 
that 

Af cos(AT + J\f)- sin (AT + AT) « -AT-AT 3 /3 . (4.1) 
Hence from (|A4|) , we expect that Eq. (|2.25 ) is well approximated by 



AT" + Af 1 - 2Af - 2A/" 3 /3 - (3e t T A Af 2 /2 - pe l T c M « . (4.2) 

where ' designates differentiation with respect to t. The relative magnitudes of the three 
source terms containing Af in the above equation as t— > — 00, will determine which one will 
dominate over the others. Assuming that the dominant behavior of these three source terms 
is 7e w * for some constants 7 and u > as t approaches —00, we can further simplify Eq. ( |4.2| ) 
and obtain the following linear non-homogeneous equation with constant coefficients: 

A7"+AT , -2AT + 7e w ' ^0 . (4.3) 

If we had exact equality in the above equation, the complementary solutions would be 
spanned by e~ 2 * and e*. A particular solution of such an equation would be of the form 
7ie wt , for some constant 71. Since AT is bounded as t— > — 00, we expect solutions for large 
negative t to be of the form 

AT sa Ce l + 71 e wt (4.4) 
for some constant C. The following two sets of circumstances may apply: 
Case I: u > 1. 

In this case, AT ~ Ce t . In order to avoid the unbounded complementary solution e~ 2t , 
we will employ the boundary condition AT'(— M) = A7"(— M). We can assign an arbitrary 
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value, 5, to both Af and to AT' at t = —M, so long as M is sufficiently large and \5\ is 
sufficiently small. Thus we can employ a standard initial value problem-solver (like Runge 
Kunta method with local error control, as can be found in any standard software) to solve 
equation (|A4j), and integrate Eq. Q2.25p from t = —M to a large positive value of t. We can 
also integrate backwards from t = — M to a negative value of t whose absolute value is even 
larger than M. 

Case II: to < 1. 

Af ~ Ce wt . In order to avoid the unbounded complementary solution e~ 2t in this case, we 
would have to employ the boundary condition AT'( — M) = uxSf(—M). 

In all the numerical experiments presented in the following illustrative examples, the domi- 
nant behavior as t— > — oo is that of the complementary solution e 4 , so that all our examples 
belong to case I. We therefore impose the initial conditions AT'(— M) = 7T( — M) = 5, for 
some large M > and small \S\ in these experiments. Provided such conditions are met, 
when we integrate equation ( |A4| ) backward in t towards — oo, numerical solutions of AT stay 



close to the value for a long interval. This will be clearly seen in Fig. 1 to Fig. 6 below. 
Since we also have an unbounded mode e~ 2t near t = — oo, Af = acts like a saddle point. 
Hence if we continue to integrate backward in time, the numerical solution will eventually 
blow up, as one can never get rid of the unbounded mode entirely in numerical calcula- 
tions, even though Theorem [I] demonstrates that solutions exist that never deviate from 
as t— > — oo. 

Should we use conditions other than 77'(— M) = Af(—M), the magnitude of the solutions 
will become large very quickly as we integrate backward in the variable t. A larger value of 
M is needed to ensure that condition ( |A5| a) is satisfied approximately in such a scenario. 
This larger domain will increase the computational cost and decrease the accuracy of the 
simulation results. 

From the proof of Theorem [l], we know that for any value of 8 and M, a solution AT 



will settle down to a value as depicted by condition (|A5| b). Hence we expect that there are 
infinitely many solutions to Eq. (|A4|) satisfying boundary conditions ( |A5|) . But the existence 
of an infinity of solutions does not imply that there are infinitely many m values that 
these solutions approach as t—>-oo. We have carried out numerous numerical experiments 
to address this question. These experiments show that solutions of Eq.([A4]) exist, for the 
identical T4, T c , and Af, for which lim t _ j ._ 0O A/* = 0, which demonstrate that AT approaches 
limits corresponding to at least three different values of m as t— >oo. We will document some 
of the more interesting numerical results below: 

A. Experiment 1 

In this experiment we have used 

M = n ^ and Ta = Tc = t 3 exp(-3t 2 ) (l - exp(-6t 2 

1 + e M r v 
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with (3 = 2grQ = 250, where we have taken — M as the value of negative numerical infinity 
with M = 13.03423. We have set Jf = Jf' = 10~ 6 at t = -M. The result, using a standard 
initial value problem-solver, is represented in Figure 1. We observe that AT is essentially 
zero from about t = —5 to t = —25 ( |A7"(— 25)|~10 -5 ), that it hovers near the unstable 
equilibrium position AT = Qtt to within better than 1 part per 10 5 at t — 10, and stabilizes 
to its final position at AT = 7n. With reference to Eq. (|3.14 ), I = and m = 7. The winding 
number for the gauge-invariant gauge field in this numerical experiment is Q = |. 



B. Experiment 2 

We employ exactly the same data here as in Experiment 1, except that in our choice of 
numerical infinity — M in this experiment, M = 13.034; and that we have set A/" = AT = 
— 10~ 6 at t — —M. The result is represented in Figure 2. We observe in this case that AT 
is again essentially zero from about t — — 5 to t — —25 ( |AT(— 25)|~10 -5 ), that it briefly 
hovers near the unstable equilibrium position AT = 6tt in the vicinity of t = 8, and stabilizes 
to its final position at AT = 5tt. With reference to Eq. (|3.14 ), I = and m = 5. The winding 



number for the gauge-invariant gauge field in this numerical experiment is Q 



5 
2' 



C. Experiment 3 



We again employ exactly the same data here as in Experiments 1 and 2, except that in our 
choice of numerical infinity — M in this case M = 30, and that we have set A* = A* = 10 -6 
at t = —M. The result is represented in Figure 3. We observe in this case that A* is again 
essentially zero from about t = —5 to t — —25 ( |A*(— 40) | < 2xl0 -10 ), that it remains 
within approximately 1% of tc in the vicinity of t = —9, and that it stabilizes to its final 
position at A/" = 9tt. With reference to Eq.( 3.14j) , t = and m — 9. The winding number 
for the gauge- invariant gauge field in this numerical experiment is Q = |. 



D. Experiment 4 

In this experiment we have used 

= ( /°° e t )2 + (l ~ T^b) 27r and T A = T c = t 3 exp(-3t 2 )(l-exp(-6t 2 )) 

with (3 = 2gro = 250, and we have taken — M as the value of negative numerical infinity, 
with M = 20. We have set AT = AT = — 10 -4 at t — —M. The result is represented in Figure 
4. We observe that A/* is essentially zero from about t = —12 to t — —30 ( |A7"(— 30)|^10 -5 ); 
that it briefly hovers near A/* = — tc from approximately t = —8 to t — —6 to within about 
.2%, and that it stabilizes to its final position at AT = 29.6874. Since tan~ 1 (27r) = 1.41297, 
and 29.6874 — 1.41297 = 28.2744 = 97r, this experiment exemplifies an £ = 1 and m = 9 
case. The winding number for the gauge-invariant gauge field in this numerical experiment 
is Q = 5.5677. 
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E. Experiment 5 



We employ exactly the same data here as in Experiment 4, but choose M = 15. In this 
experiment, we have set TV" = TV" = — 1CT 4 at t = —M. The result is represented in Figure 
5. We observe in this case that |AT(— 25)|~10 -5 , and that it stabilizes to its final position 
at AT = 17.1214. Since tan- 1 (2 7 r) = 1.41297, and 17.1214 - 1.41297 = 15.7084 = 5vr, this 
experiment exemplifies an £ = 1 and m = 5 case. The winding number for the gauge- 
invariant gauge field in this numerical experiment is Q = 3.5677. 



F. Experiment 6 

We employ exactly the same data here as in Experiments 4 and 5, but choose M = 30. 
In this experiment, we have set AT = AT = 10~ 4 at t — —M. The result is represented in 
Figure 6. We observe in this case that |A7"(— 35)|~10~ 6 , that it hovers near AT = n from 
approximately t = —17 to t = —6 to about 3% and that it stabilizes to its final position 
at AT = 35.9699. Since tan- 1 (2vr) = 1.41297, and 35.9699 - 1.41297 = 34.557 = IItt, this 
experiment exemplifies an I = 1 and m — 11 case. The winding number for the gauge- 
invariant gauge field in this numerical experiment is Q = 6.5677. 



G. Comments on Numerical Experiments 

We note that in all six experiments, we were able to calculate backward in t from t = —M 
for a long range, and still could obtain small values of AT. This shows that our choice of 
initial conditions AT'(— M) = A7(— M) has been effective and reliable. 

With reference to Theorem |l], we note that in Experiments 1, 2 and 3, the given Af leads 
to I = for the limit specified in Assumption |A2| ; similarly, in Experiments 4, 5 and 6, the 
given Af leads to t = 1. Hence, for large t, in Experiments 1,2 and 3, AT" converges to mrc 
for some integer m, and, in Experiments 4, 5 and 6, A/ converges to tan _1 (27r) + ititc for 
some integer m, in accordance with Remark 2 in the Appendix. We observe that for each 
of two values of £, m is a different integer in each experiment. Thus the numerical results 
show that there are cases in which, for the same Af, Ta, and Tc, A* can converge to at least 
three different values of m. 

A linearization analysis shows that if m is odd, then mir behaves like a stable equilibrium 
point. However, when m is even, mrr acts like a saddle point. (That is why we have a saddle 
point behavior at t = — oo, since it corresponds to m — 0, which is even.) A good analogy 
can again be made to the damped pendulum case. The even values of m correspond to an 
inverted pendulum whose equilibrium is unstable, while the odd values correspond to the 
lowest stable equilibrium point. While it is theoretically possible for a moving pendulum to 
stop exactly at the inverted position, any slight error will prevent us from observing such a 
phenomenon in numerical calculations. (This is equivalent to observing that we can never 
completely eliminate the unbounded mode in numerical calculations.) However, in principle, 
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there is no reason why AT cannot approach an unstable limit point at which m is even, as 
t^oo. In fact, Experiments 1 and 2 substantiate this claim, since the solutions AT hover 
near 6n for such an extended interval. It is not surprising to find that the precise values of 
initial conditions that bring about the transition from 5n to 7n are sensitive to the numerical 
accuracy of the computation. In Experiments 1 and 2, a change in the numerical accuracy 
of the computation can require a 1% change in the value of M to obtain similar qualitative 
results. 
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APPENDIX A: 



In Section |TTT], the discussion is based on the properties of the solutions of Eg. ( |2.25| ). 
The theorem that establishes these solutions is proven in this Appendix. In this proof, the 
function A" and N are represented as y and w respectively, and the independent variable u 
is represented as t. 

Let Ta, Tc, and w be given C 1 (continuously differentiate) functions of t in the interval 
(—00, 00), let (3 be a given positive constant, and let £ be a given integer. We assume that 
there is a positive constant a such that 

lim e {1 ~ a)t T A = 0, lim e (1 " a) % = , lim e~ at w = , (Al) 

t^ — OO t—r — CO t— > — OO 



lim e {1+a)t T A = , lim e (1+Q) % = , lim e at (w - 2£tt) = , (A2) 

t— +00 t^oo t^oo 

and 

lim e l T' A , lim e*7^ , and lim w" exist, lim e at w' = . (A3) 

(i.e. the functions w, e*7^ and e'7c are exponentially decaying at both positive and negative 
infinity.) Consider the equation 

^ + % + 2 {w cos(y + w)- sin(y + w)} 

(A4) 

+(3e l {T A (cos(y + w)-l)-T c sin(y + w)} =0. 
We have the following theorem. 

Theorem 1 Under the assumptions ( \A1\ ) to (Kh ) on w, T A and Tc stated above, there is a 
solution y to equation on the interval (—00,00). Moreover, 

(a) lim y(t) = 2nn , and (b) lim y{t) = tan -1 (2£tt) + mn , (A5) 



where n and m are integers and tan 1 : R 



— 1 . I ) ^ / 7T 7T 

2 ' 2 



Remark 1: By using the substitution Y = y — 2nir, it suffices to prove the above theorem 
when n — 0. 

Remark 2: When I = 0, then equation ( A5|b) yields lim 4 ^ 00 ?/(t) = mn for some integer m. 
When £7^0, then equation (|A"5|b) yields lim t ^ 00 ?/(t) = tan~ 1 (2£7r) + miT for some integer m. 

Remark 3: Actually the proof shows that there are an infinite number of solutions to our 
problem. 

Proof: Because of the Remark 1, we let n = 0. As t — » —00, we expect that the contribu- 
tions from the w, T A and Tc terms are negligible due to the assumptions (|A1[), so that y = 
behaves exactly like the unstable equilibrium point in the damped pendulum case. We now 
prove that this intuition is indeed what happens. 
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We rewrite the equation as 



y" + y'-2y + f(t,y) = 0, 

where 

fit, y) = 2 {(y - sin(y + w)) + w cos(y + w)} + fie* {T A (cos(y + w)-l)-T c sin(y + w)} . 

It is noted that \f(t, 0)| < o(l) e at and f y (t,y) = o(l), when t is large and negative and \y\ 
is sufficiently small. Next we convert it to a system of two equations, 

y' = z , 

z' =- z + 2y-f{t,y) . 
If, besides the conditions that / satisfies, f(t,0) also were for all t, a standard theorem in 



ordinary differential equations (theorem 4.1, p. 330, Ref. |L5| ) would enable us to conclude 
that, for any negative to with sufficiently large modulus, there is a one dimensional manifold 
S, which depends on t , in the phase space (y,y r ), such that if the initial conditions y(t ) 
and y'(t ) lie on S, then y — > as t — > — oo. In our case, we have 0)| < o(l) e a * for 
large and negative t. The same proof given in theorem 4.1 in Ref. |T5[ still works. The 
only alteration needed is in the first step, when we iterate Eq. (4.11) in that theorem, since 
f(t,0) is not zero anymore. However, because of the exponential decay of /(i, 0) for large 



negative t, Eq. (4.12) in that proof in Ref. [I5|| will still be valid, and the same proof stands. 



Fix t . We have thus constructed a solution y to equation ( |A4j ) on the interval (— oo, to) 
satisfying (|A5|a) . (Because of the one dimensional manifold S, this gives rise to Remark 3.) 



Take any M > 0. On any fixed bounded interval [to, M], 

|2 {w cos(y + w) — sin(y + w)} + fie 1 {7^(cos(?/ + w) — 1) — Tc sin(y + w)} \ 

is uniformly bounded for all y. Thus the solution y can never blow up in finite time, and 
exists for all t G (— oo, oo). 

To finish the proof of this theorem, it suffices to show that the solution has property 
(|7£5|b). For £ = 0, we expect that this is the case physically, because for large t our equation 
behaves like a damped pendulum moving under the influence of gravity only — because 
of assumptions ( |A2| ) and ( |A3| ) — hence we expect that the solution will settle down in an 
equilibrium point, i.e. y = vrm for some integer m. To prove the property (|A5| b) for arbitrary 
integer values of £, take t\ > sufficiently large so that for t > ii, w < e~ at , e*|7^| < e~ at /2, 
and e*|7c| < e~ at /2. Now for t > t\, using the Cauchy-Schwarz's inequality, 
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f t [&f- + 2 cos(y + w) + 2w sin(y + it;)] 
— y'(y" + 2w cos(y + it;) — 2 sin(y + it?)) + 2ww' cos(y + it;) 
= -(y') 2 -/3exp(t){T A (cos( 2 / + ^) - 1) - T c sin(iy + w)}y' 
+2ww' cos(y + w) 

< -{y'f + 2Pe t \y'\{\T A \ + \T C \) + 2\w\ \w'\ 

< -^qi + A(3 2 [e\\T A \ + \T c \)] 2 + 2Ce- at 

< + 8(f3 2 + C)e- at 

for some constant C > 0. In other words, if we define H = ^y- + 2 cos(y + w) +2w sin(y + 
w) + W 2 +c^ a \ then 

d -§<-^<-°- 

Hence H is decreasing in £ for t > ti, and H > —2—211^1100. Since — ^ < iJ(ti)+2+2||it;|| 00 
for t > ti, \y'\ is bounded for t E [ti,oo). From the governing equation (|A4j) , \y"\ is then 
bounded for t E [ti, oo). 

Next take the derivate of the equation (|A4|) . From this equation, we have \y"'\ is bounded 
for t E [ti, oo) due to assumptions (|A3|) . We can now conclude that both \H'\ and \H"\ are 
bounded for t E [ti, oo) by some simple calculations. 

Since H is monotone decreasing and has a lower bound, there exists a constant H such 
that lim^oo H = H . Together with the bound on \H"\, we can conclude (p. 116, Ref. HTB"! ) 
that lim^oo H' — 0. It can be checked that 



H = -W? - 0exp(t) {T A (cos(y + u?) - 1) - T c sin(y + iu)}y' 

(A7) 



+2ww' cos{y + w)- 8(/? 2 + C)e~ at . 
Taking the limit as t — > oo, we have 

lim(i/) 2 = 0, (A8) 



t— >oo 



which means y' — > as t — ► oo. Together with being bounded, this yields lim^oo y" = 
(by the same theorem in Ref. Pi). Finally, taking the limit t —>■ oo in equation (XI), we 
obtain the result 

lim (u> cos(y + w ) — sin (?/ + it;)) = . (A9) 

t — >oo 

This can be reduced to lim^oo [2£n cos y — sin y) = 0. This is equivalent to (JA5|b) . The proof 
of the theorem is now complete. 
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FIGURE CAPTIONS 



Figure 1. A solution of Eq.(|A4])) is plotted against the dimensionless variable t . The input 
into the equation is given in Subsection |1V A| . In this case, £, the winding number of the 



gauge field Af, is given by £ = 0, and the integer that characterizes the lim^oo Af, is m = 7. 
Q = | for this case. 

Figure 2. A solution of Eq. (|A4|) is plotted against the dimensionless variable t . The input 
into the equation is given in Subsection |IVB| . In this case, £, the winding number of the 



gauge field Af, is given by £ = 0, and the integer that characterizes the lim^oo Af, is m = 5. 
Q — | for this case. 

Figure 3. A solution of Eq. ( |A4| ) is plotted against the dimensionless variable t . The input 
into the equation is given in Subsection [IV C| . In this case, £, the winding number of the 



gauge field Af, is given by £ = 0, and the integer that characterizes the lim^oo Af, is m = 9. 
Q = | for this case. 



Figure 4. A solution of Eq.(A4) is plotted against the dimensionless variable t . The input 



into the equation is given in Subsection |IV D| . In this case, £, the winding number of the 



gauge field Af, is given by £ = 1, and the integer that characterizes the lim^oo Af, is m = 9. 
Q = 5.5677 for this case. 



Figure 5. A solution of Eq. ([A4|) is plotted against the dimensionless variable t . The input 
into the equation is given in Subsection [IV E . In this case, £, the winding number of the 



gauge field Af, is given by £ = 1, and the integer that characterizes the lim^oo Af, is m = 5. 
Q = 3.5677 for this case. 



Figure 6. A solution of Eq. (|A4| ) is plotted against the dimensionless variable t . The input 
into the equation is given in Subsection [IV Fj . In this case, £, the winding number of the 



gauge field Af, is given by £ = 1, and the integer that characterizes the lim^oo Af, is m — 11. 
Q = 6.5677 for this case. 
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